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Abstract 

In this article we have discovered a close relationship between the (algebraic) Bethe 
Ansatz equations of the spin s XXZ model of a finite size and the g-Sturm-Liouville prob- 
lem. We have demonstrated that solutions of the Bethe Ansatz equations give rise to the 
polynomial solutions of a second order q-diffcrcncc equation in terms of Askey- Wilson oper- 
ator. The more general form of Bethe Ansatz equations and the mathematical problems rel- 
evant to the physics of XXZ model are discussed. Furthermore, the similar correspondence 
between Bethe Ansatz equations of XXX model and the Sturm-Liouville type difference 
equation in terms of Wilson operator has also been found. 
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1 Introduction 



Among the many one-dimensional integrable quantum spin chains the XXZ spin chain is an 
important and renowned one. It corresponds to the 6- vertex model in the 2-dimensional solvable 
statistical mechanics. In particular, the XXZ Hamiltonian of spin ^ of a finite size 2N has been 
the subject of extensive studies for a long time in the physics community, and in recent years 
it has been investigated by mathematicians in the areas of mathematical physics and quantum 
algebras. The one-dimensional [/g(sZ2( C))-invariant XXZ model of spin ^ of a size 2N with 
the open (Dirichlet) boundary condition is described by the following Hamiltonian [21], [18], 

2N-1 _ -1 1-1 

41z = - E + ^^1+1 + - ^^(^? - ^In), a := (1.1) 

where are the Pauli matrices acting on the j*'* site: 

(o) 2Af 2 

The Hamiltonian Hj^j/^^ defines a linear endomorphism of ® C , whose eigenvalue problem 
has been the main concern for the physical quantities related to the system. In the context of 
quantum inverse scattering method developed by the Leningrad school in the early eighties (see, 
for example [8], [17]), the diagonalization problem of the Hamiltonian has been investigated by 
means of solutions of the following algebraic Bethe Ansatz equations, 

27V 

(sm{Xk + ^v)\ _ Y] sin(Afc + Xj + r?) sin(Afc - Xj + r?) 





^ sin(Afc -Iv)/ j4j=i ~ ^^"^'^^ ~ ~ ' 



In general, for a given positive half-integer s, the theory also provides a similar XXZ Hamilto- 
nian of spin s of a finite size 2A'^ with the following form of Bethe Ansatz equations associated 
to the corresponding Hamiltonian, 

sin(Afc + .S7?) ^ " sin(Afc + Xj + ??) sin(Afc - Xj + rj) ^ < ^ m 2) 

sin(Afc - ST]) J jj^j=i sin(Afc + Xj - rf) sin(Afc - Xj - rf)' - - ' 

For over half a century the Bethe Ansatz has been a very useful tool in obtaining important 
informations of the XXZ-model in physics literature; but yet there is no systematic study of the 
mathematical content on the Bethe Ansatz available at this time. In this paper, we establish 
an explicit connection between the relation (1.2) and g-Sturm-Liouville problem. In fact, we 
study the Bethe Ansatz equations in a more general setting than (1.2), namely the system of 
equations 

TTT sin(Afc + Sir]) ^ -pj- sin(Afc + Xj + 7?) sin(Afc - Xj + rf) ^ < < ^ m^-. 
f}^ sin(Afc - Sir]) jj^^j^^ sin(Afe + Xj - r]) sin(Afc - Xj - r])' 

where s;'s are 2N complex numbers. The mathematical problem is to obtain the solution(s) 
of the above system of nonlinear equations. The solution of (1.3) will determine the roots of 
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a polynomial which satisfies a ^'-difference relation which is a Sturm-Liouville type equation 
involving the Askey- Wilson operator, as we shall see in §4. For N = 2, the system of equations 
(1.3) is solved by the zeros of the Askey- Wilson polynomials. An interesting physical problem 
is to understand the large A'^ behavior of solutions of the system of equations (1.3). Such a 
study is a challenging problem in the area of q-Sturm-Liouville equations. The understanding 
of the solutions of (1.3) and their limiting distribution will have a profound impact on the 
physics of statistical mechanics. 

As g — > 1, it is known that the XXZ chain of spin s becomes the s/2( C)-invariant spin 
s XXX chain of a finite size L = 2N with the periodic condition. The antiferromagenetic 
spin I XXX chain of size L with periodic condition is given by the following Heisenberg XXX 
Hamiltonian: 

L 
3=1 

The above XXX spin chain is a famous integrable model with many applications to solid state 
physics and statistical mechanics. It was first proposed by Heisenberg in 1928, [11], then solved 
by Bctlic in 1931 [6]. The spectral problem of Hxxx can be reduced to the solution of the 
following Bethe Ansatz equations, 

We shall discuss the above equations subject to certain symmetry conditions imposed on the 
roots Xj's from physical considerations of the ground state of Hamiltonian Hxxx- In this 
situation, the Bethe Ansatz equations are closely connected to a Sturm-Liouville problem 
involvong a Wilson operator. The solution of this problem requires new developments in the 
theory of Wilson operators. Further mathematical study of those difference equations could 
enrich our understanding of both the mathematical and physical contents related to the Bethe 
Ansatz equations of XXX model. This will be the subject of a future work. 

An important contribution of this paper is to point out that the Bethe Ansatz equations are 
variations on nineteenth century work by Heine, Stieltjes and Hilbert. Heine studied polynomial 
solutions to a second order differential equations 

U{x)^+^x)^+r{x)y = 0. (1.5) 
ax ax 

In (1.5), n and $ are given polynomials of degrees N and N — 1. Heine proved that given a 
nonnegative integer n, there exists at most (^^"2^) choices of the polynomial r in (1.5) such 
that (1.5) has a polynomial solution. Stieltjes continued this research and showed that if we 
assume that H and $ have only real and simple zeros and their zeros interlace, then there are 
precisely i^^^'^) polynomials r which will make (1.5) have a polynomial solution of degree n. 
For references and details wc refer the interested reader to Szego's book [22] , which also treats 
Hubert's work on the location of zeros of Jacobi polynomials P^"'^^ (x) when the conditions of 
orthogonality, namely Re a > —1, Re /3 > —1, and a + f3 real, are not necessarily assumed. 
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If we denote the zeros of a polynomial solution of (1.5) by xi, . . . ,Xn then when x = Xk, for 

1 < A; < n, equation (1.5) becomes 

V 1 = ^(a^fc) ^ 

.<kn.m 2n(x,)- 

The system of algebraic equations (1.6) is a system of generalized Bethe Ansatz equations. 
Observe that although the polynomial r does not appear explicitly in the system (1.6), it 
is used implicitly to show the existence of a polynomial solution to (1.5), hence the zeros 
of the polynomial y solve (1.6). Conversely if (1.6) has a solution xi,...,x„, then we set 
y(x) = nfe=i(^ ~ ^k) and observe that Il{x)y" + ^{x)y' vanishes at the zeros of y hence there 
is a polynomial r of degree (at most) N — 2 such that (1.5) holds. This shows that the number 
of different solutions to (1.6) is the same as the number of choices of the polynomial r in 
(1.5). It is exactly this set up that is behind the modern Bethe Ansatz equations where second 
order differential equations are replaced by second orderequations in the Askey- Wilson operator 
(XXZ model) or the Wilson operator (XXX model). 

The paper is organized as follows. In Section 2, we indicate the equivalence of (1.6) to 
second order equations in the Askey- Wilson operator, which is a system of equations that 
generalize the Bethe Ansatz equations for the XXZ model. In Section 3 we provide intuitive 
explanations for closed form polynomial solutions to second order differential equations with 
polynomial coefficients of degrees 2, 1, and zero, as well as similar equations where derivatives 
are replaced by applications of g-difference and Askey- Wilson operators. This explains where 
the big Q'-Jacobi polynomials and the Askey- Wilson polynomials come from. 

Section 4 is devoted to studying the general second order equation in the Askey- Wilson 
operator T>q with general polynomial coefficients. We identify the symmetric form of such 
operator equation through an amazing simplification resulting from expanding the coefficients 
in Chebyshev polynomials of the first and second kinds. These representations lead in Section 
5 to the concept of regular singular points of the second order operator equation in an Askey- 
Wilson operator with polynomial coefficients. 

In Section 6, we consider the case corresponding to g = 1 in Section 4. The relationship of 
the Sturm-Liouville problem in terms of Wilson operator and the Bethe Ansatz equations of 
the XXX model for the ground state has been found. 

Convention. In this paper, Z, R, C will denote the ring of integers, real, complex numbers 
respectively, N = Z>o, R* = R \ {0}, C* = C \ {0}, and i = 

2 q-Sturm-Liouville Problems 

Given a function /, we set f{z) = f{x) with 

z = e'^, x = cose = {z + z~^)/2. 



Indeed z = x i: \J x'^ — \ with the branch of the square root chosen to make \/x^ — 1 ~ .x, as 
X — oo. This makes \z~^\ < 1 < \z\, and = holds if and only if x G [— 1, 1]. With this notation 
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we now introduce the following operations on f{x), 

Mix) = Rq'^'z), 



(I _ x2)-V2 

(^./)(^) = i(gl/2_g-l/2) (^./-^.-^/)(^)^ (2-1) 



The operator Vg is called the Askey- Wilson operator [9], [12]. It is important to observe that 
both Vq and Aq are invariant under q ^ q~^- 

Recall that the Chebyshev polynomials of the first and second kinds, respectively, are 

r„(cos^) = cosn^, UJcosO) = '^^^ . (2.2) 

smt* 

The Askey- Wilson operator has the properties 

„ri/2 _ „—n/2 i 
VqTnix) = _ Un-l{x), AqTn{x) = - [g"/^ + g'^/^] r„(x). 

Thus Vq reduces the degree of a polynomial by 1 while Aq preserves the degree of a polynomial. 
Furthermore the Askey- Wilson operator has the following properties, [12], 

T^.ifg) =M{Vqg) + {r^q-^g){Vqf) 

= iAqf){Vqg) + {Aqg){Vqf). 

The second line in the above equation follows from the first by interchanging / and g in the 
first line then taking the average of the two answers. 

Thus Vq,Aq are operators of the polynomial algebra C[a;] with degVq{f) = degf — 1, and 
deg^g(/) = deg/ for / G C[a;]. Furthermore these operators preserve R[.t] for a real q. 
For convenience, we shall use the following convention throughout this paper, 

^ = 6^''?, e = 2\, (hence x = cos 2 A). (2.4) 

For given functions w{x), p{x), r{x), we consider the following g-Sturm-Liouville equation of 

m, 

-^Vq Mx)Vq) /) (x) = r{x)f{x). (2.5) 
By (2.3), one can rewrite the equation (2.5) in the following form, 

Ii{x)Vlf{x) + ^{x){AqVqf){x) = r{x)f{x), (2.6) 



where the functions 11, $ are defined by 



II{X) = -^Aqp{x), ^{X) = -^Vqp{x). (2.7) 

w{x) ^ ^ ' w{x) 
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The form (2.5) is the symmetric form of (2.6), which can be seen from the formula of integration 
by parts in [7] which will be stated later as (5.3). We shall show in §4 how to construct w and 
p from n and Now the relationship 

(A^./)(x) ^ 



(g-l)(ge2ie-l)(e2i«-g) 



Km 



2(73/2 



i(l - g)2 sin 0(ge2i« - l)(e2i« - q) 



x{(e'^^ - g)v + (^e^'' - " + " l)}/(^), 



shows that a root xq = cos2Ao of the polynomial f{x), f{xo) = 0, necessarily satisfies the 
following equation, 

{(e^^^° - q)iU{xo) - $(xo)sinr/sin2Ao)v}/(a^o) 

+|(ge4iAo _ i)(n(xo) + $(xo)sin?7sin2Ao)r/g-2}/(xo) = 0, 

or equivalently, 

f Vq^f \ ^^^^ ^ - sin(2Ao + r?) (n(xo) + ^(xp) sin r? sin 2Ao) g) 



\r]q-2f J sin(2Ao - ri){U{xo) - $(a;o)sin77sin2Ao 

For a polynomial f{x) of degree n with distinct simple roots xi, . . . ,Xn, one writes 

n n 

/(x) = 7 n (a; - xj) = 7 n (cos 2A - COS 2X/) ,7/0. 
It is straight forward to see that 

n 

^ 3=1 



_ g-K^+-^fc--^j)) _ g-i{'?+-^fc+Aj)^gi(??+Afc-Aj) _ g-i(j?+AA:-Aj)^| 

n 

(-1)^7 n sin(Afe + Aj + vi) sin(Afc - Aj + 77); 



n 

= {-'^Tl n sin(Afc + \j - rj) sm{\k - \j - rf). 

Now observe that (2.8) indicates that the roots xi,...Xn of f{x) satisfy the system of 
equations, 

-sin(2Afc + r?)(n(xfc) + ^(xfc) sinr/ sin 2Afc) 
sin(2Ajt -ri){U(xk) - ^>(xfc) sin?7sin2Afc) 

_ sin(Afc + Xj + T]) sin(Afc - Xj + r/) ^ < < „ 
^ sin(Afe + Aj - 77) sin(Afe - Xj-rj)' 
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In other words we arrive at the system of nonhnear equations 



= n ^|n(A, + A, + ,)sin(A,-A,+,) ^^^^ 



n(xfc) + sin T] sin 2Afc 

n(xfe) - $(a;fc) sin 77 sin 2Afc j4,i=i ~ '^^ ^^^''^'^ ~ ~ 

which we call the Bethe Ansatz equations associated with Il{x),^{x). For the interest of 
applications to physical problems, we shall only consider those g-Sturm-Liouville problem where 
the coefficients Il{x),^{x),r{x) of (2.6) are polynomials in x with the degrees 

degn =l + deg$ =2 + degr =N>2. (2.10) 

3 Second Order Equations 

We first introduce some notations from [2], [9]. The g-shifted factorials are defined by 

n 

{a;q)o = l, {a; q)n =]]_{!- aq''~^), n = 1, . . . , or 00. (3.1) 

fe=i 

Furthermore 

Formula (3.2) clearly holds when a is a nonnegative integer but is used to define the general 
g-shifted factorial when a is not necessarily an integer. 
Askey and Wilson [4] introduced the polynomials, 

^„(cos6';a) := (ae'^,ae"'^;g)„ = (ae'^; g)„(ae~'^; g)„ , (3.3) 

as a basis for the space of polynomials. This is the most suitable basis for use here. A different 
basis appeared in the g-exponential functions in [15] where it was used to provide a g-analogue 
of the expansion of a plane wave in spherical harmonics. Clearly 

(j)^(x; a) = (-2a)"g"("-^)/2x'* + lower order terms. (3.4) 

It is straightforward to see that 

lim cf)n{x; a) = (1 - ae'^)"(l - ae"'^)" = (1 - 2ax + a'^)''. 

We may use (3.2) to define the more general functions {(j)oi{x;a)} by (3.3) when a is not 
necessarily an integer. It readily follows that 

Vq^a{x; a) = (l)a-i{x; aq^l'^), 

Aq(t>a{x;a) (x;agV2)[l_ag-V2(l+^a)^ + a2^a-l]. 

The second formula in (3.5) holds when a = provided that we interpret the product defining 
(f)a in (3.3) as in (3.2). Furthermore we have 

2AqMx; a) = (1 + q-")ct>a{x; aq'/^) + (1 - g-")(l + a\'"'-')<j>a-i{x- aq^'^). (3.6) 
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We shall use TTjix) to denote a generic polynomial in x of degree j. Consider a differential 
equation 

Tr2{x)y"{x) + 7ri(x)y'(x) + Xy{x) = 0, (3.7) 

where A is a constant. We seek a polynomial solution to (3.7) of degree n. We know that one of 
the coefficients in tti or 7r2 in not zero, hence there is no loss of generality in choosing it equal 
to 1. Thus TTi and 7r2 contain four free parameters. The scaling x ^ ax + b oi the independent 

variable absorbs two of the four parameters. The eigenvalue parameter A is then uniquely 
determined by equating coefficients of a;" in (3.7) since y has degree n. This reduces (3.7), 
in general, to a Jacobi differential equation whose polynomial solution, in general, is a Jacobi 
polynomial, see [22]. Solutions also include special and limiting cases of Jacobi polynomials 
including the Bessel polynomials and the plynomial x^. 
Next let us consider the same problem for the operator 

= M^IM. (3.8) 

X — qx 

Consider the operator equation 

Tr2{x)D''gy{x) + Tri{x)Dgy{x) + Xy{x) = 0. (3.9) 

Here under the same assumptions on tti and tt2 one easily finds out that the only scaling allowed 
on X is X ^ ax, hence one of the coefficients in tti and it2 is chosen as 1 and the remaining 
coefficients of tti and 7r2 constitute three free parameters. The general polynomial solution to 
(3.9) is the big qf- Jacobi polynomial which contains three free parameters, [16], [9], [2]. 
We next consider the Askey- Wilson case 

'K2{x)Vly{x) + Tri{x)AqVqy{x) + Xy{x) = 0. (3.10) 

In this case one can not perform any scaling on x, so apart from assuming that one of the 
coefficients in tti and 1^2 is unity, we have four free parameter, namely the remaining coefficients 
in 1T2 and the coefficients in tti. This is the case of the Askey- Wilson polynomials where 

7r2(x) = [2(1 + CT4).x2 - ((71 + a3)x - 1 + £72 - £74] , 

7ri(x) = 2 [2(^4 - l)x + (71 - CT3] /(I - g), 

where aj is the jth elementary symmetric function of parameters ai , 02 , 03 , 04 . In order to 
solve (3.10) for y we expand y in the Askey- Wilson basis {(t)n{x;a\)} and find 

A = -4g(l - q-^){l - a^q^-^){l - q)'^. (3.12) 

With the above choice for A the polynomial solution to (3.3) is unique and is given by an 
Askey- Wilson polynomial of degree n, sec for example [7]. 

At this stage one wonders whether replacing 7r2, tti and A by vr^'+i, vr^, and iTk-i in equations 
(3.7), (3.9) and (3.10) lead to more general orthogonal polynomials. Grunbaum and Haine 
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[10] proved that the only orthogonal polynomial solutions to (3.10) after the replacements 
(7r2,7ri,A) — > (TTfc+i, TTfe, TTfc-i) are the Askey- Wilson polynomials or special and limiting cases 
of them. Ismail [13] showed that the same conclusion holds without assuming orthogonality. 
This generalizes earlier work of Hahn, and Bochner who proved that the g-Jacobi polynomials, 
and the Jacobi polynomials are the only polynomial solutions to (3.9) and (3.7), respectively 
with the above replacements. 



4 Multiparameter Cases and Bethe Ansatz Equations for the 
XXZ ]V[odel 

For 2A'' complex numbers aj,! < j < 2N, we denote 

a = (ai, . . . ,a2Ar), 

and aj the j-th elementary symmetric function of a^'s for < j < 2N, with ao := 1. We define 
the weight function w{x,a), 

( (ei^^,e-^^^gf)o. 

- sin M n|^,(a,ei^ a,e-i^; qU ' ^ 
which can also be written in the following form, 

_ 2ie-f^(e'^^gf e-'^^gf ),, _ -2ie^(gf e^^^ e"^^^; gf ),, 
''^'''''^ ^,^^l(%ei^a,e-i^;g)oo ^,^fl(a,■ei^ a.-e-^^; g)oo ' ^ ' ^ 

With w{x) = w{x,a),p{x) = w{x,q^/'^a) in (2.5), we shall consider the following equations, 

—^Vq{{w{x, q^'^a)Vq)y){x) = r{x)y{x), r{x) e C[x] , deg r = N - 2. (4.3) 
With the notation 

n(x; a) = / ^. Agw{x, q^^^a) , $(x; a) = / ^. Vqw{x, q^^^a), (4.4) 
w{x,a) w{x,a) 

the equation in (4.3) becomes 

U{z; a)Vly + $(z; a)AqVqy = r{x)y. (4.5) 

Note that for N = 2, r(x) = r G R and aj,q G R with |aj| < 1,0 < g < 1, the weight 
function wg{x) is a positive function on (—1, 1), and the solutions of the equation (4.3) are the 
Askey- Wilson polynomials which are orthogonal in L^((— 1, 1); (1 — x^)~^/^), see e.g. [7]. 

Theorem 4.1 The functions Il(x;a),^{x;a) are polynomials of x of degree N,N — 1 respec- 
tively, and have the following explicit forms, 

N-l 

U{x;a) = -q-m{(^-l)^aM+^{-iyi<Ji + a2N-i)TN-iix)}, (4.6) 

1=0 

$(x;a) = _ E(-l)'(^'-^2iv-0^iv-;-i(^)- (4.7) 

^ ^ 1=0 
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Conversely , for given polynomials 11 and $ of degrees N and N — 1 respectively, there is a 
unique 2N -element set {aj : 1 < i < 2N} such that n(a;) = n(x, a) and $(a;) = $(a;;a) and 
(4-4) holds. 

Proof: By the definition of w(x,d), also the forms in (4.2) when applying riq,r]q-i, Il{x;a) has 
the following expression, 



iq^ Sin ^ ni=i(«,ei', a.e"'''; g)oo 



which implies 



e 2 (g 2 e , e 



e 2 (e'"", g 2 e 



a) = \q 4, sm 



2A'' iM. 2N 

' ^ n(i-v^)-r^n(i- 



1 - e 



-Ar/4 



2Af 2Af 

j=i j=i 



i=i 

-i9 



2N 



i=0 



JV-l 



-gVK-l)^'^^ + E + a2N-i)TN-i(.x)}. 



By the same method, we have 



NB 



^{x; a) = 



2q— sin ^ 
(gi/2 _g-i/2)sine 



-i^'g 27V live 2Af 

n (1 - + 11(1 - a,e-'^) 



2Af 



2Ar 



(gV2-g-i/2)sine 



(4.8) 



2N 



(gl/2 _ g-1/2) Z^^ ^ 



sin^ 



2rr 



N-l 



— Im 13 (-l)'(c^/ - (T2N-i)Un-i-i{x). 



(^l/2_g-l/2) . 

To see the converse statement, given H and $ we expand them in Chebyshev polynomials of 
the first and second kinds respectively, then define ao by cro = 1 and aN by (— times 
the constant term in the expansion (4.6) in terms of Chebyshev polynomials. Then define the 
remaining cr's through finding ai it a2N-i from the coeflficients in 11 and $ in (4.6) and (4.7). □ 
It is important to note that Theorem 4.1 gives a constructive way of identifying the param- 
eters ai, . . . a2N and w{x; a) from the functional equation (4.5). 
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Theorem 4.2 Let ai = q-^' = e^^*"^' for 1 < I < 2N. The Bethe Ansatz equations (2.9) 
associated with the polynomials n(x; a), a) have the form (1.3). 



Proof: By Theorem 4.1, (indeed in its proof), we have 
n(a;; a) = 



-N 

-q 4 



2 

-iV 

$(xfc;a)sin?7sin2Afc = 
hence 



2N 2N 

g-2iiVA, a,e'^^^) + e^'^^^ - a^e-'''^^' 



2N 2N 



2N 

i=i 

2N 

n(xfe;a) - $(xfc;a)sinr?sin2Afc = -gT^e-^'^^'^ JJ(1 - a^e^^^*). 

By substituting aj = q~^^ in the above formula, the result of this theorem follows from (2.9). 
□ 

If \q\ > 1 replace q by 1/p, use the invariance of Vq and Aq under q q~^ to rederive 
Theorem 4.2 with q replaced by 1/p. This covers the case \q\ > 1. The cases \q\ = 1 and in 
particular the cases when is a root of unity do not seem to be amenable to the techniques 
developed here. 

The Bethe Ansatz equations (2.9) describe the relations of roots Xj,l < j < n, of a 
polynomial f{x) of degree n in the g-Sturm-Liouville problem (4.3). It is important to note 
that for given 11 and $ in (2.6) with degll = N and deg$ = N — 1, there are 2N complex 
numbers, ai, a2, . . . , a2N, such that I^{x) = n(x; a) and $(x) = ^{x; a). For a positive integer 
N and all aj = q~^ with s G ^Z>o, the equation (2.9), or (1.3) become (1.2) which is the Bethe 
Ansatz equations for the spin s XXZ model of a size 2N with the open (Dirichlet) boundary 
condition [21]. This correspondence shows that identifying the spectrum of XXZ spin chain is 
related to spectral problems of a g-Sturm-Liouville equation. 

For N = 2 , we have Il{x;a) = 7r2(x) and ^{x;a) = 7ri(x) in (3.11). In the special case 

ai = — a2 = e~^*, 03 = —04 = e~'^*''"''\ 

the polynomial solution to (4.3) is a g-ultraspherical polynomial of degree n. In this case (1.3) 
becomes 

sin(2Afc + s) sin(2Afc + s + rj) _ -pj- sin(Afc + Aj - ij) sin(Afc - Xj - rj) ^ < < „ 
sin(2Afc - s) sin(2Afe - s - r/) f |_ sin(Afe + Aj + r?) sin(Afe - Aj + r/) ' ~ ~ 
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For N = 2 with Sj < for all j, y{x) is an Askey- Wilson polynomial of degree n. In fact, one 
has the following result. 

Theorem 4.3 Let N = 2 and r] = > 0. Then for all n the system (1.3) has a unique 

solution provided that sj < 0, 1 < j < 4. Furthermore all the X's are in (0, 7r/2). 

Proof. Let y be a polynomial of degree n with zeros cos(2Aj), 1 < j < n. Wc know that (1.3) 
implies the validity of (3.10) for x = cos(2Aj). Here tti and tt2 are as in (3.10). With A in (3.9) 
chosen as in (3.12) the left-hand side of (3.10) is a polynomial in x of degree n — 1 and vanishes 
at n points. Hence (3.10) must hold for all x, and y must be an Askey- Wilson polynomial of 
degree n. Since the Askey- Wilson polynomials are orthogonal on [—1, 1], all their zeros are in 
(-1,1). □ 

For large n the distribution of the cos(2Aj)'s follows is an arcsine distribution. This follows 
from general theory of orthogonal polynomials since in this case j\ \nw{x){l - x'^y'^/'^dx is 
finite [19], [20]. Note that the proof of Theorem 4.3 used the orthogonality of the Askey- Wilson 
polynomials in an essential way. 

Using arguments that go back to Heine, one can proved that for given polynomials H and 
$ of degree N and N — 1 respectively, there exist at most {^^^2"^) choices of r{x) such that 
(2.6) has a polynomial solution of degree n. This shows that ioT N > 2 the solution to (1.3), 
without additional assumptions, may not be unique. 

For ttj = q~^i the monic Askey- Wilson polynomials satisfy the recurrence relation 

XPn{x) = Pn+l{x) + (cos 2?? + 2A„ + 2C„)p„(x) -I- AAn-lCnPn-l{x), (4.9) 

where 

_ s\n{{n-l-a[)ri)Y{j^2^M{n-sx-Sj)r,) 
" ~ sin({2n-l-a'^)ri) sin((2n-o-^ )r;) ' 

(4.10) 

_ sm(nr7) f]2<j<fc<4sin((n-l~Sj-Sfc)»7) 
~ sin((2n-l-cr'Jr;) sm((2n-2-cr'Jr;) ' 

and a'j is the jth elementary symmetric function of the Sfc's, (see (3.1.4)-(3.1.5) in [16]). From 
the form of and C„ in (4.10) it is clear that the positivity condition A„_iC„ > may not 
hold for all n. 

Theorem 4.4 Assume that ij and all the Sj 's are real and that AqCi > 0. Let 

M = Max{n : Ak-iCk > 0, for 1 < A; < n}. 
Then (1.3) has a unique solution for < n < M and the X's are all real and simple. 

Proof. We have a family of M -|- 1 (which may be = oo) orthogonal polynomials and their zeros 
are all real and simple. It is known that in this case (4.5) has only one polynomial solution [5]. 
□ 
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Remark: It is important to emphasize that when the orthogonahty condition ^„_iC„ > is 
violated, then the Askey- Wilson polynomials continue to satisfy (4.5), hence their zeros solve 
the Bethe Ansatz equations but wc can no longer guarantee the reality or the simplicity of the 
zeros of the Askey- Wislon polynomials. □ 



5 Singularities and Expansions 

In this section we consider different regimes. We assume 0<g<lsoir7<0 only in (5.1)-(5.4) 
below. Recall the inner product associated with the Chebyshev weight (1 — on (—1, 1), 

namely 



</,5>:= j'j{x)g{x)-0=. (5.1) 



For < g < 1, one observes that the definition (2.1) requires f{z) to be defined for |gf='=^/^2;| = 1 
as well as for I z| = 1. In particular Vq is well-defined on H1/2, where 

■= {/ : /((^ + l/z)/2) is analytic for < \z\ < g"""}. (5.2) 

Brown, Evans, and Ismail [7] proved the integration by parts formula 



<Vqf,g> = 



- < /, Vl^Vq{g{x) (1 - x')-^/'') >, 



(5.3) 



for f,g£ Hi/2- They also proved that if wlx) > on (—1, 1), and y G Hi, p & Hi then the 
eigenvalues of the g-Sturm-Liouville problem, 

1 



w[x 



-Vg {p{x)Vgy{x)) = \y{x), (5.4) 



are real and the eigenfunctions corresponding to distinct eigenvalues are orthogonal with respect 
to w. This implies the following theorem. 

Theorem 5.1 Assume that Q < q < 1, \aj\ < g^/^, 1 < i < 4, and w{x\a) > on (—1,1) 
where a = (01,02,02,04). Then the eigenvalues of 



^-^Pg (w{x,q^/'^a)Vqy{x)j = Xy{x), 



w{x, 

are real and simple, and the eigenfunctions are mutually orthogonal on (—1,1) with respect to 
the weight function w{x; a). 

Proof. The condition |oj| < q^^"^, 1 < j < 4 ensures that p{x) = w{x;q^^'^a) is in Hi and the 
theorem follows from [7]. □ 

We now discuss the concept of a regular singular point of the equation (4.3), which is 
the general form of a second order equation in the Askey- Wilson operator with polynomial 
coefficients as we say in Theorem 4.1 for a general q and A'^. 
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First recall that the concept of singularities of differential equations is related to the analytic 
properties of the solutions in a neighborhood of the singularities. We have no geometric way to 
describe the corresponding situation for equations like (4.3). In the present set up the analogue 
of a function analytic in a neighborhood of a point (a + a~^)/2 is a function which has a 
convergent series expansion of the form X^^o '^4'n{x', a). We have no other characterization of 
these g-analytic functions. 

In the case of second order differential equations the singularities are the zeros of the 
coefficient of y" and the regular singular points (^j 's are precisely those points where one can 
construct a series solution to the differential equation of the form J2^=o'^niz — (^j in a 
neighborhood of z = (j. A closer examination of the equation (4.3) in the form (4.5) reveals 
that one can formally expand a solution y as X^^o Uni'nix; a), substitute the series expansion 
in (4.5) and recursively compute the coefficients y„ provided that a is not among the 2N 
parameters {d, . . . , C2iv}, where C,j = {aj + aJ^)/2, and possibly the points and oo, which we 
do not know how to handle. This indicates that what plays the role of singular points of (4.5) 
are Ci) • • • C2N in addition to a: = 0, a: = 00. Around the singular point x = (^j, we shall expand 
functions in the set {(pn+a{x', aj)} , where (l)a{x;a) is as in (3.3). We let 

00 

yi^) = J2y'^ (pn+a{x; ai), (5.5) 
n=0 

and observe that r{x)(f)n+a{x', oi) is a linear combination of {(f)rn+a{x', oi) : n < m < n+N — 2}. 
Furthermore we note that (3.2), (3.3), and (4.4) imply 



r^-^Pg (w{x; q^^'^d)T>q(pn+aix; 01 
x: a) \ 



w{x; a) 

2ai{q — l)w{x; a) 
1 - 

2ai(g — l)w{x; a) 



Vg \ w{x;q^^'^a)(l)n+a-i{x;q^^'^ai 

Vg (w;(x;aig"+'*-V2,a2,...,a2iv)) 



(l-g"+")^;(x;aig"+" «2, ■ ■ ■ , a2jv) . n+a-i „ „ x 

— — — — <P{x, aiq ^ , 02, • • • , a2N) 

2ai{q — l)w{x; a) 

■ $(a;; aigr''+""\ 02, ... , a2N) (pn+a-i{x; ai). 



2ai{q-l) 

The substitution of the expansion (5.5) for y in (4.3), reduces the left-hand side of (4.3) to 

2a ( -1) 02, • • • , a2N) 4>n+a-i{x; ai)yn. (5.6) 

n=0 ^ 

Note that the smallest subscript of a in r{x)y{x) in (4.5) is a. On the other hand (5.6) 
implies that appears on the left-hand side of (4.3). Thus the coefficient of (pa-iix',ai) 

must be zero. To determine this coefficient we set 

AT-l 

$(x;g"~^ai,a2,...,a2jv) = ^ dj{q"') (t)j{x;aiq"'~^), (5.7) 

j=0 



14 



and after making use of (f)n{{ci + cl ^)/2; a) = 6nfi we find 

do(g") = Hiaiq"-' + ar'g'-")/2; (?"-^ai, 02, • • • , a2Ar). 

Thus the vanishing of the coefficient of <pa-i{x;ai) on the left-hand side of (4.3) imphes the 
vanishing of (1 — q°')do{q'^), that is 

(1 - (/")^.((ai(7"-^ + ar'g'-")/2; q"-'ai,a2, a2Jv) = 0. (5.8) 

Theorem 5.2 Assume \aj\ < 1, for all j. Then the only soluUon(s) of (5.8) are given by 
qa ^ 1^ 0^ qa ^ q/^aiaj), j = 2,...,2N. 

Proof. Prom (5.8) it is clear that = 1 is a solution. With x = (aig"~^ + ai^q^~°')/2 as in 
(5.8) we find e'^ = aiq°'~^, or ai^q^~°'. In the former case, 2isin0 = aig"~^ — ai^q^~°', hence 
(5.8) and (4.8) imply 

which gives the result. On the other hand if e*^ = q^~°'/ai, then we reach the same solutions 
via (4.8). □ 

Ismail and Stanton [14] used two bases in addition to {(f)n{x; a)} for polynomial expansions. 
Their bses are 

Pnicose) := {l + e'''){-q'-^e^'';qX_,e-'^^, (5.9) 
0„(cos^) := iq'/^e^',q'/^e'^';q'/X- (5-10) 



They satisfy 



V,pn{x) = 2g(i-")/2l-^p„_i(a;), (5.11) 



1 - , 

-DgM^) = -2gV4i^^^_^(a,) (5.12) 

One can also seek solutions of second order operator equations by expanding solutions in the 
above polynomial basis. This will be the subject of future work. 

Atakishiyev and Suslov [5] studied certain expansions of solutions of a very special non- 
homogenous equation corresponding to (3.10), that is iV = 2, with special value of A. They 
did not however investigate any concept of singularities, nor they have observed the general 
structure of expanding the polynomial coefficients in Chebyshev polynomials. Their Wronskian 
is inadequate because according to their definition, the Wronskian of two polynomials is not a 
polynomial. Askey and Wilson wrote down the functional equation (3.10) and identified it as 
the equation satisfied by the Askey- Wilson polynomials. Of course when one works only with 
the case N = 2, as in [5], one is bound to miss the complications and the elegance of the case 
of general N. 
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6 Bethe Ansatz Equations for the XXX Model and Wilson Op- 
erators 



In this section we consider the problems arisen from the discussion in the previous sections 
when q tends to one, the corresponding Sturm-Liouvillc problem connecting with the Bethe 
Ansatz equations of the Heisenberg XXX spin chain. We write the variable z in previous 
sections in the form, 

z = e'^ = q-'y, 

and again the parameters, aj = q~^L It is known that as q tends to 1, the transformation f(z) 
and the operators riq±i,Vq,Aq become f{y),r]±,W, A, respectively, where 



fiy) 

{Wf){x) 



= f{x) with X = y ; 

= ^^M+f -'n-f){x) , 

= i(77+/ + r/_/)(x). 



The above divided difference operator W is called the Wilson operator [24]. We have the 
relation 

W{fg) = {Wf){Ag) + {Af){Wg). 
Analogous to the g-Sturm-Liouville problem (2.5), we consider the following difference equation 

^ -W {p{x)Wf) {x) = r{x)f{x) , 



w{x) 

which is equivalent to the Sturm-Liouville problem in the form, 

Uix)W''f{x) + ^x){AWf){x) = r{x)f{x), (6.1) 
where 11, $ are the functions defined by 

U{x) = -^Ap{x), ^x) = ^Wpix). (6.2) 

For our purpose with reason which will be clearer later on, we seek polynomial solutions f{x) to 
(6.1), but n(x), ^{x),r{x) are rational functions with the same degree constraints as in (2.10). 
From the relations 

AW fix) = j^^iiv - ^)vl -{y + + i}/(a^); 
w'm = ^^){{y - h)vl + {y + - 2y}f{x), 

it follows that if /(xq) = for xq = y^, then 

^ -(yo + ^)(n(xo)-<i'(xo)7/oi) 
[ef i (yo - 4)(n(xo) + $(xo)yoi) ' ^' ' 
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For a degree n polynomial f{x) with roots Xj for 1 < j <n, and let x = j/^, xj = that is 



Then 



r]\f{xk) = 7 n (yfc ~ =^ + =^ 



By (6.3), ?/j's satisfy the following system of equations, 
-(yfc + |)(n(xfc)-i/fci$(xfc)) 



(2/fe-|)(n(xjt)+yfei^(xfe)) 



f=i (yfe - Vj - i)(yfe + Vj - 1) 



or equivalently, 



n (^^-^^• + ^)(^^' + ^^- + ^), l<k<n. 



n(a;fc) + yfei$(xfc) ^-^^ j=i {Vk - Vj - i)(2/fe + Vj - i) 



(6.4) 



We now consider the equation (6.1) which arises from (4.3) by letting g — 1. In the notation 
of the g-gamma function [2] , [9] 



r,(y) 



we make the identification 

1 



{q";Q)oo 

^ ^q{-Sj-iy)^q{-Sj+iy) 
(l-g)2+2«.(g;g)2. 



, aj = e . 



Note the limiting property lim^^i Tq{y) = r(y), for a proof see Appendix I in [1]. The limiting 
weight function (4.1) has the following expression: 



w{x, a)dx 



The fact limg_ 
ing to g = 1, 



(1 - q^f{q^]q^)l, VS=i ^gjsi - T^y)Tq{-si + iy) logg ^ ^ 
(1 - g)'^+'S.^^(g;g)4A^(-iiVy)r . (iTVy) sin ^ 



dy. 



sin y 
'0 sin(iV0/2) 



2/N identifies the following weight function as the one correspond- 



w{x, s) 



n^r(-s, + iy) 



r(iiVy) 

Consequently the function p{x) in the Sturm-Liouville problem (6.1) is given by 

p{x) = w{x,s-'^) 1 := (i,...,i). 



(6.5) 



The corresponding functions in in (6.2) will be denoted by 
n(x; a) 



—^—Aw{x,s- \), ^{x; a) = — — - 
w[x,s) 2 w[x^s) 



Ww{x,s ) 

2 



17 



Theorem 6.1 For a given s = {si, . . . , S2n) with an even N, denote i^j the j-th elementary 
symmetric function of Si's for < j < 2N, {<^o := 1). Then Il{x;a),x^{x;d) are the polyno- 
mials of X of degree at most N with following expressions, 

Uix:a) = (-l)f {z^ + Ef=on-iy(^<?2,+i - ?2,+2)x^-^'-H, 
x^x;a) = (-l)f {Ef=-o^(-iy (1,2, - + h2N}. 

The roots, Xk = y"^, k = 1, . . . ,n, of a degree n polynomial solution f(x) of the Sturm-Liouville 
problem (6.1) satisfy the following Bethe Ansatz type relation, 

yk^fivk + sj}^ -j^ {m-y, + ^){m + yj + -^) ^^^^ 

yk + ^ ,ii vk - sii j^tj=i ~ yj ~ ^>^y^ + y^ ~ ^> 

Proof. It is easy to see that 

. nj=ir(-^,-ij/ + i)r(-.s,+iy) 

2^ - T{-iNy + ^)TC.Ny-^) ' 

, . 1, ngir(-gj- iZ/)r(-5, + iy + l) 

ri-wix,s 1 = — Tf Tf . 

' ^ ' 2^ r(-iAry-f)r(iiVy+f) 

The functional equation of the Gamma function, r{z + I) = zr{z), establishes the explicit 
representations for 11 and 

y j=i j=i 

N -. N-1 

= (-l)^{^(-iy,2,x^-^- + I E (-iy^2,-MX^-^-^} 

j=0 ^ j=0 

N-1 



j=Q 

$(x;a) = i-^{(y + l)n(y-is,)-(y-^)n(2/ + ^)}; 

y ^ j=i j=i 

N-l . N 

= (-i)^{ E {-iy^\2,+ix^-^-' + 1 Ei-^y^2,x^-'-'} 

j=0 s=0 

N-1 . , 

= (-1)^{ E (-1)^(2^2, - ?2,+l)x^-^-^ + ^<^2NX-'}. 
j=0 

Hence we obtain the expressions of n(a;; a), a). By the first expressions for Il(x; a), a) 

in terms of y in the above right hand sides, one can easily derive the fohowing identities: 

(~i)f ■ 2^ f— Tit ■ 2^ 

U{x;a)-^x;a)yi = ^ — —{y-^) Uiy+isj), U{x;a)+^x;a)yi = ^ — —iy+^) Uiy-isj), 

y ^ j=i y ^ j=i 
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hence 

11(3:; g) - <^{x;a)yi ^ Vk - ^ y + sji 
U{x; a) + ^>(x; a)yi + 5 /=i V - Sji ' 

By the above relation and (6.4), we obtain (6.6). □ 

Remark. For the Bethe Ansatz equations (6.6) with N odd, one can reduce the problem to 
the above theorem for some even N' by adding certain zero- value Sj's. By the similar method, 
one enables to apply the above theorem to Bethe Ansatz problem of the following type with 
si, . . . , SM e C and M G N, 

l[yji±ll= fr im-y, + i){m + y, + -^) 

□ 

For a positive half-integer s, it is known that the spin s XXX model of an even size L with 
the periodic boundary condition has the following Bethe Ansatz equations: 

(^)^= fl A.GC, fe = l,...,/. (6.8) 

At. — SI . , , , Ajt. — Ao — 1 

In the antifcrromagenetic case, the ground state is on the sector I = i^; conjecturally there is 
the unique real solution of the Bethe Ansatz equations. Note that {—Xj^j^i is a solution of 
(6.8) whenever {Aj}^-^j^ is a solution. Hence for the ground state, the roots AjS are expected 
to be real and invariant under the sign-change ( up to permutation of the indices j). Having 
this ground state conjecture in mind, we now consider a general problem of (6.8) for I = ^ 
(mod 2) with roots invariant under the change of sign, \j ^ — Aj, i.e., A^s with the following 
form: 

L = 4M + 2, {A,}^.=i = {0,±yi,...,±y.„}; (6.9) 
L = AM, {XjYj^i = {±2/1, . . . , ±y„} . (6.10) 

In this situation, one can link the Bethe Ansatz equations (6.8) of XXX model to the Sturm- 
Liouville problem previously discussed in this section. In the case (6.9), the relation (6.8) 

becomes 

{yk-\){yk--^){yk + s\Y^^^ _ " (y, - t/^ + i)(i/, + j/^ + i) ^ 

ivk + '^{yk + i){yk - siY^+^ - (% -yu- i)(%- - + i) ' ^ " ^' " " " ' 

in which case, we set the iV, Sj's in Theorem 6.1 as follows: 

L = 4M-|-2, s^l: N = 2M + 2, si = 0, S2 = -1, S3 = ■ ■ ■ = S2N = s, 

^2N = 0, ,j = C^r')s^ - Cf-,')s^-\ < i < 2Ar - 1; 
s = l: N = 2M+1, si = 0, 53 = • • • = S2Ar = 1, 

^2iv = 0, 9 = ('^^"') 0<j<2Ar-l. 

(6.11) 
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In the case (6.10), one has 



n 



{vj -yk + i)iyj + + 
ivj -yk- '^){yj -yk + i) 



, Vke C\{0}, k = l,...,n, (6.12) 



j=l,¥k 



in which case, N, Sj in Theorem 6.1 are given by 



L = 4M, N = 2M, si 



S2N 



'2iV' 



for < j < 2N. 



(6.13) 



For both situations, the relations are of the form (6.7) for suitable s;'s. Note that the non-zero 
condition of y^s in (6.9), (6.10) for the corresponding polynomial solution f{x) of (6.1) requires 
one further constraint, namely /(O) ^ 0. The ground state of antiferromagenetic spin s XXX 
model of size L is governed by the real root solution of the above equations for n = ^. 

We now consider the case s = ^ and discuss some mathematical problems relevant to the 
physics of the system. The Hamiltonian is given by Hxxx in (1-4). We will discuss the Bethe 
Ansatz equations (6.12) with s = ^ in (6.11) (6.13). To illustrate the mathematical content 
connecting to the Bethe Ansatz equations for the Hamiltonian Hxxx, we consider the case 
L = 2, 4. For L = 2 in (6.9), by 



|r(iy)r 



y sinh iry ' 
the weight function is expressed by 



r(^ + iy) 



4-IT 



{4y^ + 1) cosh Try ' 



w(x) = 



r(iy)r(l+iy)r(^+iy)2 



(27r)3y 



(x + 1)^ sinh27ry 



r(2i2/) 

By (6.11), the corresponding Sturm-Liouville problem (6.1) is governed by the equation 

(8x2 + 6x + l)W^f{x) + {Ax + l)AWf{x) = Xf{x) 

with A = 4n(2n — 1) except n = 1, in which case A = 4. For L = 4 in (6.9), by (6.13) we have 
the following weight function and the Sturm-Liouville problem. 



w(x) 



r(^ + iyr 



r(2iy) 



Air^y sinh ivy 
{x + l)"^ cosh^Try 



-x{32x^ - 16x - 6)W^f{x) + (48^2 + 8x- l)AWf{x) = (Ax + /x)/(x) 

with A = — 16n(2n — 1) except n = 1, in which case A = 48. For n = 1, we have = 4 and 
f{x)=x + j. Note that due to the parameters in the above weight functions, the Sturm- 
Liouville problems we encounter here does not make the Wilson polynomials orthogonal, as 
in [4]. However, our empirical computation for the case L = 2 has indicated that the /(x)'s 
are real polynomials with /(O) ^ 0. For the general size L, the weight function of the Sturm- 
Liouville problem connecting to XXX model in the context of discussions in this section has 



20 



the following expressions: 



L = 4M + 2 



w{x) 




L = 



AM. 



w{x) 



(x + cosh'^ ny 



Nevertheless a sound mathematical treatment of their corresponding Sturm-Liouville equations 
appears still lacking now, and it remains a difficult task to obtain substantial knowledge in- 
corporating the physical applications associated to XXX model. However The quest of a such 
mathematical theory would be a necessary one in order to understand the essential features of 
the Bethe Ansatz equations for the Heisenberg XXX spin chain. 

7 Conclusions and Perspectives 

In this paper we have provided a one to one correspondence between polynomial solutions 
to Sturm-Liouville type equation involving the Askey- Wilson operator and solutions of the 
algebraic Bethe Ansatz equations of XXZ and XXX models. We have also started a preliminary 
study of the mathematical problems arising from the physics of the Bethe Ansatz. In doing 
so we have reduced the nonlinear problem of solving the Bethe Ansatz equations to the linear 
problem of finding polynomial solutions to certain linear second order equations in terms of 
the Askey- Wilson or Wilson operator. Thus the physics of XXZ model has indeed raised some 
mathematical questions which demand the need for a systematic mathematical development 
in the theory of the corresponding q'-difference equations. The mathematical solution of those 
questions will lead to a better understanding of some interesting physical quantities of XXZ 
and XXX models. We intend to continue this study in future and partial results obtained so 
far are indeed very promising. 
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